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NON-LOCAL EFFECTS BY HOMOGENIZATION 
OR 3D-1D DIMENSION REDUCTION 
IN ELASTIC MATERIALS REINFORCED BY STIFF FIBERS 

ROBERTO PARONI* AND ALI SILlt 


Abstract. We first consider an elastic thin heterogeneous cylinder of radius 
of order e: the interior of the cylinder is occupied by a stiff material (fiber) 
that is surrounded by a soft material (matrix). By assuming that the elasticity 
tensor of the fiber does not scale with e and that of the matrix scales with e 2 , 
we prove that the one dimensional model is a nonlocal system. 

We then consider a reference configuration domain filled out by periodi¬ 
cally distributed rods similar to those described above. We prove that the 
homogenized model is a second order nonlocal problem. 

In particular, we show that the homogenization problem is directly con¬ 
nected to the 3D-1D dimensional reduction problem. 
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1. Introduction 

Quite often nature combines two, or more, materials with complementary prop¬ 
erties to generate a structural material whose performance and functionality su¬ 
persede those of monolithic materials, [25] ■ Natural materials usually present a 
complex hierarchical structure, with characteristic dimensions spanning from the 
nanoscale to the macroscale, that is quite difficult to replicate in fabricated materi¬ 
als, [26]. Nowadays, engineers and material scientists design and fabricate novel 
materials, inspired by natural materials, with a complex, but not hierarchical, 
structure. Two materials with completely different properties are appropriately 
arranged to create advanced functional materials that are lightweight, with re¬ 
markable strength, stiffness, and toughness, ®mm • The material properties of 
a two component material can be tuned by appropriately choosing the individual 
components, their morphology, their size, and arrangement, [51. The understand¬ 
ing of the mechanical response of these advanced functional materials is an issue of 
paramount importance for industrial applications. 

In the present paper, we study structures made up by a stiff and a soft material; 
the latter usually provides energy dissipation, toughness, ductility, and makes the 
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structure lighter. This combination is quite common in engineering and biological 
composites. 

We first derive, within the framework of linear elasticity, the elastic problem 
governing the motion of a rod composed of a stiff and a soft material. This is 
achieved by considering a sequence of problems posed on cylindrical reference con¬ 
figurations of diameters that scale with a parameter e, and by taking a variational 
limit as e approaches zero. The interior of the cylindrical regions is occupied by the 
stiff material (fiber), while the surrounding part by a soft material (matrix), whose 
modelling is achieved by scaling the elasticity tensor by e 2 . In this way the ratio 
between the components of the elasticity tensor of the fiber and those of the matrix 
is 1/e 2 . No assumption on the material symmetry of the body is made and the 
material could be inhomogeneous also within each of the stiff and the soft regions. 

We show that the limit problem may be written equivalently as two independent 
systems: the first posed on the matrix, while the second takes into account the 
clastic energy of only the fiber and the loads applied to the entire body. Mathe¬ 
matically speaking, the problem posed on the matrix is local, while the problem 
posed on the fiber is non-local since also the loads applied outside of the fiber 
region enter into the problem. For any distribution of loads the problem is still 
non-local, since in order to get the global information about the limit displacement 
the problems posed on the matrix and on the fibers have to be both solved. Leaving 
technicalities aside, the problem on the fiber determines a Bernoulli-Navier type of 
displacement, while the problem on the matrix determines the deviation from the 
Bernoulli-Navier type of displacement that takes place within the matrix region. 
The need of a displacement correction within the matrix region could be explained 
by the fact that the matrix is much more deformable than the fiber. 

As a first step, to find the variational limit, we deduced a priori bounds on the 
displacements. These are easily obtained within the fiber region by means of Korn 
inequality, while a more intricate argument is needed within the matrix region. 
Indeed, since there is a loss of uniform ellipticity within the matrix region, only 
through a judicious application of several Poincare and Korn inequalities, that also 
exploit the a priori bounds obtained within the fiber region, we are able to prove a 
priori bounds within the matrix region. 

We also address a related homogenization problem. We consider a reference 
configuration domain Q filled out by periodically distributed rods similar to those 
considered in the rod problem mentioned above. To completely fill , we take 
rods with a square base of size e. The body is therefore made of two regions: one 
occupied by a stiff material (fibers), while the second is filled by a soft material 
(matrix); the ratio between the two elasticity components is still 1 /e 2 as in the rod 
problem. We pass to the limit by using two-scale convergence, this introduces a 
slow variable x and a fast variable y. Both variables appear in the limit problem: 
the x variable is used to describe the macroscopic behaviour while the y variable the 
local, or microscopic, behaviour near each fiber. The limit problem could be recast 
in terms of only the macroscopic variable x, by taking advantage of the linearity of 
the problem; we refrain of doing it for the sake of brevity, see 0 . One of the main 
features of the limit problem is that locally each fiber behaves like the rod problem 
previously considered. As a consequence the homogenized problem is a second order 
problem: we recover a second-order material as a variational limit of a first order 
material, see also jH |20]. We finally emphasize that the techniques used connect 
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the homogenization problem directly to the 3D-1D dimensional reduction problem. 
In particular, the non-locality already found in the rod problem still appears in the 
homogenization problem, for exactly the same reasons. 

The literature on homogenization and on dimension reduction for the diffusion 
or elasticity problems is huge, see for instance m 0 0 1 nu m um na na on m . 
We here describe only the works that are more related to ours. 

In [4] the homogenization of a periodic medium composed by stiff fibers sur¬ 
rounded by a soft matrix, similar to that considered in this paper, is studied. Both 
the fiber and the matrix are taken to be homogeneous and isotropic. Also, the 
elasticity tensor of the matrix scales with e 2 , like in our case, but the scaling of the 
elasticity tensor of the fibers is proportional to 1/e 2 . The limit problem is non-local 
but it does not depend on the extensional energy of the fibers. 

The same homogenization problem with the same scaling as in |4], but in a more 
general context including anisotropic and inhomogeneous materials, was considered 
in [9]. In addition to the non-local phenomenon arising in the limit and already 
proved in [4], it was shown in |9] that the anisotropy of the fibers occurs in the 
limit problem in the same way as in |15j and m leading to nonstandard terms. 
The phenomenon due to anisotropy occurs also in the present setting but to show 
it we would have to rewrite the problem in a reduced form. For the sake of brevity 
we do not give such formulation. 

We have chosen to scale the elasticity tensor of the matrix by e 2 and to keep 
the elasticity tensor of the fibers fixed, i.e., it does not scale with e, among all the 
possible scalings, in order to have a non vanishing extensional energy of the fibers 
in the limit problem, differently from what happens in Si- 

Also in [5] a similar homogenization problem, still for isotropic and homogeneous 
matrix and fibers, has been considered. The diameter of the fibers is allowed to 
scale with size different from e, the elasticity tensor of the matrix is kept fixed, 
while that of the fiber may take several scalings. By taking the scaling of the latter 
equal to 1 /e 2 and by multiplying the full problem by e 2 our scaling is recovered but 
only with a load that scales at least as e 2 : that is, a sequence of loads equivalent 
to a zero load. 

In both papers HE], the fibers and the matrix are taken to be homogeneous 
and isotropic, while we do not make this assumption in the present paper. It is well 
known, see 0 HU EDS [24] that when the material under consideration presents high 
heterogeneities the nature of the limit problem is quite involved and a particular 
nonlocal phenomenon may appear in the limit. 

In the next section we fix the notation and we formulate the problems. In 
Subsection 12.11 we define the rod problem while in Subsection 12.21 we present the 
homogenization problem. In Section [3] we state our main results and we make some 
remarks. Section [4] and Section [5] are devoted to the proofs of the main results for 
the rod and homogenization problems, respectively. 


2. Notation and setting of the problems 

As usual, and unless otherwise stated, Greek indices a,/3, 7 ,..., take values in 
the set {1,2} while Latin indices i, j,k ,..., in the set {1,2,3}. Also, the summation 
convention for repeated indices is adopted throughout the paper. 
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Often, for convenience, a 3 x 3 symmetric matrix E is written as 


(1) 


f E a p E a 3 \ 

\ E a3 E 33 J ’ 


where, with an abuse of notation, we denote by E a p the 2x2 matrix whose com¬ 
ponents are E a p, and similarly E a3 denotes the vector whose two components are 
E a 3 . The product of two such matrices E , K is defined by 


( 2 ) 


E ■ K := EijKij 


( E a p 

E' a 3 \ 

( K a p 

K a 3 \ 

V E a3 

E 33 j 

' V K a3 

k 33 j 


A generic point x £ M 3 is denoted by x = (x',x 3 ), where x' = ( X\,X 2 ) £ R 2 . 
The notation x R indicates the point obtained from ( x 3 ,X 2 ) by a rotation of tt/2 
around the ^ 3 -axis, i.e., x R := (—X 2 ,Xi). The gradient with respect to x and x' 
are denoted by V and V / , respectively, and the strain of a differentiable function u 
is denoted by 

Eu := symVw := -(V« + Vu T ). 


Throughout the paper <9,; denotes the partial derivative with respect to a 

In the homogenization problem we shall have, besides the variable x, also a 
variable y. We denote by 



for any differentiable function <p. Also, 


( 3 ) 


{E u <P) a p ■■= ^(dlvp 


d V 8<Pa), 


( E yx < p )«3 :=2»3 


d 3 ^ a ), 


for every differentiable vectorial field <p. 


2.1. The rod problem. Let u>,u> C M 2 be two open, bounded, simply connected 

O O 

sets with Lipschitz boundaries, duj and dui, respectively, such that ui CC u>. 

O 

We assume the x a -axes to be centered in the center of mass of to so to have that 
(4) / x a dx' = 0. 

J CO 

Let £ > 0, I := (0,£), and 

Sl:=wxl, F:=u>xl , M:=Q\F, 

where F denotes the closure of F. Physically, we can think of as the reference 
configuration of the rod, F as the region occupied by the fiber (strong material), 
and M as the region occupied by the matrix (soft material). If x € S2 we have that 
x' £ uj and x 3 G /. 

For every e > 0, let R e £ R 3x3 be the diagonal matrix whose entries are e, e, 
and 1, i.e., R e = diag(e,e, 1). The scaled gradient H e u and the scaled strain E e u 
are defined by 

H e u := (R^-^uiR*)- 1 , E e u := (R E )~ 1 Eu(R e )~ 1 = sym H £ u. 

In components we have 

(E u)a.p y^/3(FI ri) a 3 d 3 u ai 
£ £ 

(H e u) 3 p = \d p u 3 , 


( 5 ) 


(■ H s u ) 33 = d 3 u 3 . 
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Let C be a fourth-order symmetric tensor field defined on 17. We assume that C 
fulfills the following assumptions: for all i,j,k,l, for all 3 x 3 symmetric matrices 
E, and for almost all x G 17 


{ u) (x) — ^- r jikl{x) — ^-'klij(x')i 

b) c ijkl € L°°(n), 

c ) there exists m > 0, C ijkl(x)E k iEij > m EijEij. 

We assume the rod to be clamped at both ends of the cylinder 17, we thus set 

H\ d {n) := {U G H\n-R 3 ) : u( s',0) = u(x',£) = 0} , 

and subjected to body forces / G L 2 (f7;K. 3 ). In the notation above the two d’s in 
H\ d are used to recall that on both ends of the beam we impose Dirichlet boundary 
conditions. 

By u £ G H dd (£l) we denote the solution of 


( 7 ) 


uS e Hdd.(P) 

[ (Xf + £ 2 XM)C(x)E e u e ■ E E ip dx = [ f(x) ■ ip(x) dx Vtp G H dd (VL). 

J n J n 


By the previous assumptions on the tensor C and the body force /, problem 0 is 
well-posed by the Lax-Milgram Theorem. 

Problem 0 is a rescaled elasticity problem, as explained in Remark l‘2.II below, 
of a body whose elasticity tensor is equal to C on F and equal to e 2 C on M. Thus, 
the elasticity of the fiber, i.e., the elasticity of the material occupying the region F, 
is 1/e 2 larger than the elasticity of the matrix (material in M). We point out that 
our assumptions allow both the fiber and the matrix to be inhomogeneous and to 
be fully anisotropic. 

One of the aims of the present paper is to deduce the variational limit of 0 as 
e goes to zero. 


Remark 2.1. Problem 0 is the variational formulation in the fixed domain 17 of 
an elasticity problem posed in the variable thin domain f7 e , see |10j . Indeed, let 

I7 £ := eui x /, F e := euj x /, M e := fl e \ F e . 

Let r £ : 17 — > I7 e be defined by r e (x) := {ex\,£X2,x^). Then R e = V?’ e . For any 
v : 17 —>• M 3 we define v : I7 e —>• R 3 by 

v := (R £ )~ 1 v o (r E ) -1 , 

in components 

Va(£Xi,£X2,X3) = -V a (x), V 3 (£X 1 : £X2,X 3 ) = V 3 (x). 

£ 

We note that 

Vt) = Vu o (r £ )~ 1 (R £ )~ 1 = H £ v o (r 6 )" 1 , Ev = E £ v o (r 6 )" 1 , 

With this notation, problem 0 rewrites as 
( 8 ) 

r u £ &H\ d (i\) 

] [ (\F e + £ 2 XM e )C £ (x)Eu £ ■ E(pdx = [ f £ (x) ■ <p(x)dx, \Fp £ H\ d ({ 7 e ), 

l. J n e J £i e 
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where 

C £ -Co^)" 1 , f e := R e fo(r e )~ x . 


From the previous definition, the components of the load have two different scalings: 
( 9 ) fa =£f a o(r e )~\ fi = f 3 o(r*)-\ 


2.2. The homogenization problem. The reference configuration of the elastic 
body is now assumed to be the parallelepiped 11 = w x (0,f) = w x J, where £ > 0 
and u) is the square u> = (—£,£) 2 . We consider a sequence e approaching zero such 
that 

w = U w h ere w* = eY + e(*i, if) and Y = (— -, -) 2 . 
ieh 

Let 12* and I e be defined by 

^ := u>l x (0, £) 


and 


Is ■= {* = € z 2 , w* C w}. 

Let D be the disk centered at the origin with radius 0 < r < Set 


D\ := eD + 


The set F e of fibers contained in 12 and its complementary set in 12, the matrix M e: 
are defined by: 


( 10 ) 


F e := |J FI where F* := D\ x (0,f) 
iei e 

< M e := Ct\F e = (J Mg, 

ieie 

12 = \J 12* = F e U M e . 

, »e/ £ 


The microscopic variable will be denoted by y = (j/i, j/2) (notice that there is no 
fast variable in the vertical direction) and, as before, y R = (—j/2, J/i)- 
We now assume that the elasticity tensor C is such that: 

( a) C ijki(x,y) = Cjiki (x, y) = C k uj(x,y), 
b) Cijki £ L°°(12; C#(F)), 

c) there exists m > 0, C ijki{x,y)E k iEij > m E i: jE ij: 


for almost all (x,y) £12 x F, for all i,j,k,l , and for all symmetric 3x3 matrices 
E, and where C#(F) denotes the space of continuous F-periodic functions on R 2 . 

As for the rod problem, we assume the body to be made up from a set of strong 
fibers surrounded by a soft matrix, with ratio between the elasticity tensor of the 
soft part and the strong part being equal to 1/e 2 . Therefore, the variational problem 
modeling the equilibrium of the heterogeneous clastic medium at the microscopic 
level may be written as: 


(12) 




uS £ Hddityi ' 

[ ( XF e + £ 2 XM e )C(x : —)Eu e ■ 

Jn £ 



£ fa(x, ^-) • <p a (x) 
x f 

+ f3{x, —) ' <f 3 (x) dx, 
£ 


. Vy> £ 
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Note that compared to (|8|). only the shape of the domain fl has changed. Note 
also that we have assumed, in m, the loading f(x,y) to depend on both the 
macroscopic variable x and the microscopic variable y : and that the scaling of the 
loads coincides with that of (|8jl. see ©• 

By assuming fi £ L 2 (fi; C#(Y)) : we deduce the existence and the uniqueness of 
the solution u e , for each fixed e > 0, by means of the Lax-Milgram Theorem. 

To derive the limit problem, we will use compactness properties related to the 
two-scale convergence, see [ 1 ] and El- Recall that a sequence u £ £ L 2 (fl) (weakly) 
two-scale converges to a function u £ L 2 ( x F), written u £ —^ u , if 

/ u e (x)ip(x, —) dx —> / u(x,y) ip(x,y) dxdy V ip £ L 2 (fl; C#(F)). 

J n £ JQ.xY 

A sequence u e £ L 2 (VL) is said to strongly two-scale converge to u £ L 2 (fl x Y), 


written u £ 


it, if 


u e (x)v £ (x) dx —>• / u(x,y)v(x,y) dxdy 
JnxY 


iQxY 

for any v £ £ L 2 (Vl) that weakly two-scale converges to v £ L 2 (fi x Y). 

We recall a few properties of the two-scale convergence that will be useful in the 
analysis below: 

• any bounded sequence in L 2 (VL) admits a two-scale converging subsequence; 

• if u £ is bounded in iJ 1 (n) then there exist u £ H 1 (J7) and u 1 £ L 2 (fl; H#(Y)) 
such that, up to a subsequence, 

,i. 


V X U + Vylt 


u £ —^ u, W 

if u e is bounded in L 2 (fl) and if eVu E is bounded in T 2 (f2; 
exists u 1 £ L 2 (n-H^(Y)) such that, up to a subsequence, 


then there 


eS7u e 


V y u L 


if u E £ L 2 (VL) weakly two-scale converges to u £ L 2 (fi x Y) and 


lim ||u e 
£->0 


Il 2 (0) 




L 2 (QxY) 


then 


u, see 


3. Main results 

The main results of the paper are stated in the present section. 

3.1. The rod problem. We first define a few spaces that will be useful in the 
results stated below. The space of Bernoulli-Navier type of displacements satisfying 
the prescribed boundary conditions is denoted by 

BN dd (F) := {u £ H^F) : (Eu) ia = 0 a.e. in F} 

= {n £ H^(F) : 3£ Q £ H 2 ( 0, £), 3^ £ ^(0, £) s.t. 

u a = £ a , U'i = £.3 - x a £' a , a.e. in F}. 

We note that the elements of BJ\Tdd{F) may be naturally extended to elements of 
BNddity- The set of functions orthogonal to the bi-dimensional rigid displacements 
is denoted by 

7 ZV^iF) := {(wi,w 2 ,0) : w a £ T 2 (/; H 1 (w)), / x\w 2 — x 2 w\ dx ' = 0} 

J UJ 
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while the set of twisting and cross-sectional warping displacements is 

Tldd(F) := {(-x 2 '0,x 1 ’d,v 3 ) : with i? e Hq(0,2) and v 3 G L 2 (/; H^u))}. 
where 

H L&) ■= {v G H\u) : -f dx' = 0}, 

J UJ 

For brevity, we set 

U := BAf d d(F), V:=H dd (F), W:=KV^{F) : 

and 

Z := {z : Zi £ L 2 (/; i? 1 (w)), Zi = 0 a.e. in F}. 

The following theorem is our main result for the rod problem. 

Theorem 3.1. There exists (u,v,w, z) £ U x V x W x Z such that the sequence 
of solutions u e of problem 0 fulfills the following convergences 

—> u a in U 1 (r 2 ), 

U% -X u 3 + z 3 in L 2 (I;H 1 ( w)), 

pe-.e., —i. ( (-^ W ) a P {Ev) a3 \ . t2(q\ 


E £ u £ xm 


(Ez'jctp 2 d a z s \ . 7 - 2 /r\\ 

ia„ 2s o )» m L < n> ' 

The limit ( u,v,w,z ) is i/ie unique solution of the problem 
(13) 

(u,v,w, z) € U x V x W x Z, 


C(x) 


( Ew) a0 (Ev) a3 
(Ev) a3 (Eu) 33 


{Ew) a/3 {Ev) a 3 
(-E«)e*3 (^u)33 


C(ar) 


(Ez'jcf} 2^ a ^3 
\d a z 3 0 


Xf dx + 

Xm dx = 

fa(x)u a (x 3) + f 3 (x)(u 3 (x) +Z 3 (x))dx, 


\d a z 3 0 


V(it, v,w,z) £U x V x W x Z. 


Remark 3.2. Some remarks are in order. 

• The space WxVxWxZ equipped with the norm 

II (' U,V,W,Z ) || 2 = 71(11 ( Eu ) 33 ||^,2(Q) + || ( Ev ) a3 ||^2(Q) + || (Ew) a p ||^2(Q) 
o:/3 

+ || (Ez)af 3 || 2 2(Q) + II d a Z 3 || 2 2(Q)) 

is a Hilbert space (for an analogous proof in the homogenization framework 
one can refer to [5], Lemma 4.5). Therefore, the well-posedness of the 
limit problem (fT51) is easily obtained from the Lax-Milgram Theorem since 
/ £ (L 2 (H )) 3 and C is strongly coercive, 0. 
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• The limit displacement can be described, in words, as follows: it is a 
Bernoulli-Navier displacement with the axial component augmented over 
the matrix region. Said differently, Z 3 takes into account the deviation 
from a Bernoulli-Navier type of displacement in the matrix region. The 
need of this correction over the matrix region could be explained by the 
fact that the matrix is much more deformable than the fiber, because of 
the different stiffnesses. 

• The system posed on the fiber F coincides with a problem obtained in 
m this is not surprising since the sequence u e has ( E e u e )ij bounded in 
L 2 (F) : as in [15]. In particular, it was shown in [16], and also in [9] in the 
homogenization framework, that the displacements ( v , w) may be expressed 
in terms of the Bernoulli-Navier displacement u. For brevity, we will not 
reproduce here the proofs of these statements. 

• Choosing successively (u,v,w) = (0,0,0) and then 5 = 0, we see that 
problem m may be written equivalently as two systems: the first one 
posed on the matrix M with unknown z and the second essentially posed 
on the fiber F with unknown (u, v, w). Indeed, the second problem considers 
the elastic energy of the fiber and the loads applied to the fiber and the 
matrix. Thanks to the Bernoulli-Navier structure of the displacement u one 
can redefine the loads and rewrite the second problem over the domain F 
only; the redefined loads will depend on the original loads that were acting 
also on the matrix. 

The well-posedness of these two problems may be established indepen¬ 
dently, but these problems present a non-local effect: the displacements in 
the fiber are influenced also by the loads applied to the matrix. 

Also, to describe the limit displacements, even in the region M, the solu¬ 
tion of both systems are needed, as, for instance, stated by the convergence 
u% —> U 3 + z 3 in L 2 (I; F[ 1 (co)), note that Z 3 = 0 in F, see Theorem ld.ll 

The non-local effect mentioned in the previous remark is further studied in the 
following theorem, which also gives the limit of the strains defined on the variable 
domain. 


Theorem 3 . 3 . Let (u,v,w,z) be the solution of (fldl) . Let 


E , = ( (Ew) a p ( Ev ) a 3 


E m :— 


_ f (Ez'jap 

\d a Z3 0 


{Ev ) a3 (Eu ) 3 3 J ’ 

Then the sequence of solutions u e of (jSJ| is such that 

j- o Eu e dx' —> + Ef dx ', -j- ^ Eu e dx' -A j- ° E, 

J ecu Juj J e(uj\u)) Juj\oj 

in L 2 (/;R 3x3 ), 

j- eu E a dx 1 —> u a {x3) in L 2 {I), 

J £UJ 

j- u\ dx' —> U(x 3) := j- (u3 + Z3) dx' in L 2 (I). 

J £OJ J UJ 

U may be written as 

U(x3) = | « 3 dx + mo(x 3 )j- q / 3 dx' + m 0 o(x 3 ), 

J UJ J uj\uj 


dx' 
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with mo and moo given by 


vno{x z ) = -r-j f z 3 dx', m 00 (x 3 ) = A / z 3 ° dx', 

M Joj\& M J u \s, 

where z° and z 00 are respectively the solutions of the following problems 
z° £ L°°(I; Z°) 


(14) 


/ \ 

f UJ\0. 


C(x) 


(Ez°) a p \d a z% \ ( (Ez) a p \d a z 3 


I a z o 

2 u <* z 3 


0 


kd a z 3 


0 


dx = 


/ z 3 (x')) dx' a.e. in I 

J uj\lu 


. \/z € Z° := {z : Zi € ff 1 (w \ u>), Zi = 0 on dco}, 


(15) 


' z 00 £ Z, 
[ C(a:) 

JM 


Vz £ Z. 


{Ez 00 ) af i ±d a z$° \ ( ( Ez) a p \d a z 3 \ 

\d a z\f 0 J ' V \d a z 3 0 ) dl 



f 3 (x) dx')z 3 {x)dx 


Remark 3.4. 


• The last convergence in Theorem 3.3 emphasises the non-local effect. This 
can be easily seen in the particular case in which the loading term f 3 de¬ 
pends only on x 3 . In this case, we have z 3 ° = 0 and hence 



dx' U{x 3 ) = + u 3 dx + mo(x 3 )+ f 3 dx' in L 2 (I). 
J UJ J U)\oJ 


Thus the axial displacement in the limit problem is the sum of two terms: 

— j- u 3 dx' which comes by solving a problem over the fiber but by taking 

J UJ 

into account also the loads over the matrix; 

— mo{x 3 )j- f 3 dx' which is deduced by solving a displacement problem 

J uj\uj 

over the matrix region. Note that from (fl4l) it follows that mo{x 3 ) > 0 
for a.e. x 3 £ I. 

The displacement 

U(x 3 ) = + u 3 dx' + m 0 (x 3 )+ ^ f 3 dx' + m 00 (x 3 ), 

J UJ J uj\uj 


has been split in the sum of three terms: the first is the average of the 
Bernoulli-Navier displacement, the second takes into account the correction 
due to uniform loads over the matrix region and the last the presence of 
variable loads with zero average. 

• All the previous results remain true when the cylinder fis assumed to be 
clamped on only one of its ends; the limit system is the same as m with 
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the spaces U and V defined as before but with the space H\ d (VL) replaced 
by 

:= {u G H\n;R 3 ) : u{x'J) = 0} , 

while the spaces W and Z remain unchanged, since the boundary conditions 
do not enter in their definition. 

3.2. The homogenization problem. The spaces we will use to describe the ho¬ 
mogenized problem are the following (below, h stands for homogenization): 


U h :={(u u u} a ,^) 

(16) 

V h :=L 2 (fl;<(D)), 

where 

H l m (D) := {ip G H\D) : / ydy = 0}, 

JD 

and 

W h := {(wi,w 2 , 0) : w a G L 2 ( Q; H 1 (D)), f yiw 2 - y 2 W\ dy = 0, a.e. x G fl}. 

JD 

The following theorem, analogous to Theorem 13.11 is our main result for the 
homogenization problem. 

Theorem 3.5. There exist G U h , v 3 G V h , and w G W h such that the 

sequence of solutions u e of problem m fulfills the following convergences 

eu E a ->• u a , in H 1 (Q), eVu e a V'u a + V y u\, 


: u a G H\Q) n£ 2 (w;ff 2 (/)),tf G L 2 (w; 

< € L 2 (n- H^(Y)) n L 2 (fl; H^D)) n L 2 (co x D; H 3 (I)), 
and u 3 G L 2 (n;H#(Y)),u 3 (-,y) G L 2 (lv, Hq(I)) for a.e. y G Y, 
such that 

u 3 (x,y) = -y a d 3 u a (x) + f u 3 (x,y)dy, for a.e. (x, y) G x D 

JD 

Ua(x,y) = -ypdpu a (x) +d(x)y% for a.e. (x,y) G O x D}, 


u 3 u 3 , 

Moreover, setting v a := ya^> we have 


(F ,.e w 1 ( E y w) a p (E yx v) a3 \ ( A 

(Eu ) XFe ( ( E y * v ) a3 (Eu ) 33 ) 


e 


^ f (Eu + EVu 1 )^ (.E"*ii) Q3 


(E« x u) 


o3 


e(Eu e )xM t 


Xy\d(v)■ 
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The limit (m, u*, d, v 3 , w) is the unique solution of the problem 
(17) 


D 


c 


(E*w) a p {E yx v) a z 

(■ E«*v) a3 (Eu) 33 


(EVw) aP (Ey*v) a 3 

(^ x «)a3 ^3 


dydx 





Y\D 


c 


{Eu + Eyu l ) a p (Ey*u) a3 \ 
{Ey*u) a3 0 ) 

( (Eu + E y u l ) a p ( Ev*u) a3 
' V {E yx u) a3 0 


dydx 


[ f u 

J n Jy 


x, y) ■ u a (x) + f 3 {x, y) ■ u 3 (x, y) dydx, 


V(Ui, Ua,i9, v 3 , w) € U h x V h x W h , 


where v a := yff d. 

Remarks similar to those made in Remark 13.21 still hold. 

The limit problem, even if it is not immediately apparent, is a second order 
problem since the term (Eu) 33 involves second derivatives of the function u a , as it 
can be easily deduced from the admissible set m . Thus we recover a second-order 
material as a limit of a first order material, see also [H 120 . 

One can highlight the nonlocal effect in the homogenized equation through the 
analogous of Theorem 13.31 as follows. 

Theorem 3.6. Let d,v 3 ,w) be the solution of 1171) . Let 

( (Ey W ) a p ( E« x v) a3 \ „ _ ( (.Eu + Ey u l ) a p (Ey*u) a3 \ 

f - (EV*v ) a 3 (Eu )33 )' m - V (E yx u ) a 3 0 J ■ 

Then the sequence of solutions u e of (1121) is such that 

Eu e XF e / Ef dy in L 2 (Q; R 3 * 3 ), Eu £ xm c / E m dy in L 2 (Ll;R 3x3 ), 
Jd Jy\d 


eu s a — 1 u a in i7 1 (n), 
u% —*• U(x) := j- u 3 dy in L 2 (Q). 

U may be written as 

U(x) = j- u 3 dy + m 0 (x)4- f 3 dy + m 0 o(x), 
Jd Jy\d 


where mo and moo are given by 
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being z u and respectively the solutions of the following problems 




( 18 ) 


IY\D 


C(x,y) 


(E y Z°) a p \ ( (E v z) a p §» 


\dlz% o 


dy = 


2 d a Z 3 0 

/ z 3 (y)dy for a.e. x £ fl, 

JY\D 


(19) 


„ Vz £ Zg := {z : Zi € H 1 (Y ), Zi = 0 in D}, 

C(x,2/)( 

\ 2 U OL 

= {f 3 {x,y)-f f 3 (x,y) dy)z 3 (x,y))dxdy 

Jqx(y\d) Jy\d 


lnx(Y\D) 


( Ey Z ° 0 ) af) 19^3°° 

2^3™ 0 


(Eyz) a0 \d v a z 3 \ 

\d v a z 3 0 ) dxdy 


VzeL\n ; z£). 


Similar remarks to those made in Remark 13.41 apply to the present contest after 
replacing the variable x' by y. 


4. The rod problem: proofs of the main results 

We start by stating an a priori estimate and then work up to it. 

Theorem 4.1. Let u e £ H dd (fl) be the solution of problem (|7)l. Then 

(20) sup (\\E e u e \\ L 2 ( F) + \\£H e u e \\ L 2 {F ) + IK||hi(f)) < +oo, 

£ 

(21) sup (\\£H S U £ \\ L 2 {M ) + ||(«l,wi)||fl-i(M) + ll«!ll i3(J . H i (u \5))) < +°°> 
in particular 

(22) sup (||c74 £ M e || L 2(n) + ||(«i,< +oo. 

£ 

To prove the theorem it will be convenient to set 

E e (v) := f (xf + £ 2 XM)<C(x)E e v ■ E e vdx. 

Jn 

We now state and prove several short lemmas. 

Lemma 4.2. For every v £ H dd (Li), we have 

Ee(v) > c\\E £ v\\ 2 L 2( F ), E\ s (v ) > c||i;||^-i( F ). 

Proof. Using the strong ellipticity of C, see © c), we have 

T’ £ (v) > rn\\E e v\\ 2 L 2^ > m\\Ev\\ 2 L 2^ F ) > c||n|| F i( F ), 

where the last inequality follows from Korn’s inequality. □ 
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Lemma 4.3. Let R,R + ,R C R 3 be three open regions with Lipschitz boundary 
such that R + n R~ has a strictly positive two dimensional measure and the closure 
of R + U R~ coincides with the closure of R, R. Then , there exists a constant c > 0 
such that for every v £ 74 1 (i?; R 3 ) we have that 

\\Vv\\ L 2 {R+) < c(\\Ev\\ L 2 ( R +' ) + ||Vu|| i 2 ( i J-)). 

Proof. Assume that such a constant c does not exist. Then for every integer n 
there exists a function v n £ H 1 (R-,M. 3 ) such that 

(23) ||Vu"|| i 2 (i j+ ) = 1, \\Ev n \\ LHR+) + IIWIU^*-) < i. 

Thus, ||Vi» n ||i 2 ( fl ) < 1 and the sequence v n dehned by 


:=v n --f v n dx, 
Jr 


is bounded in 


for some v £ H 1 (R; 


so that up to a subsequence we have 
v n ^v in H l (R] R 3 ), 


But, from (E3l) we infer that Vu = 0 on R~ and that 
Ev = 0 on R + . Thus v = k on R~, with k £ R 3 . As a consequence we first deduce 
that the trace of v on R + 0 R~ is k and then that v is equal to k also on R + , 

Ev = 0 on R + . Thus v = k on R. Since 7- v n dx = 0 we deduce that also 


since 


Jr 


+ v dx = 0 and hence that k — 0, i.e., v = 0 on R. By applying the standard 
Jr. 

Korn’s inequality on R + , 

ll^"'l|pf 1 (R+) < c(\\Ev n \\ L 2 {R+) + ||u" ||l 2 (r+)) 

we see that ||u n ||i/i(R+) approaches zero since the right-hand side converges to zero. 
But this contradicts (ESI). □ 


Lemma 4.4. There exists a constant c > 0 such that for every v £ we 

have 

\\H £ v\\l 2 (m) < c(\\E e v\\ L 2 ( M ) + \\H e v\\ L 2 ( F )). 

Proof. We first extend v by zero, without renaming it, to w x 1. Then we 
have that v £ 4f 1 (w x R) and that, with the notation introduced in Remark m 
v £ 74 1 (ew x R). We now partition ew x I with slices of thickness e, that is we 
write 

euj x R = [J S'*, S* := eoj x {ei, e{1 + 1)), 
ieh 

where I e is the set of the corresponding integer values of i. We also set 

F' e :=EU,x(£i,E(i +1 )), Mi~Si\F2, 

and 

S := ui x (0,1), F := ui x (0,1), M~S\F. 

Let v l : S —X R 3 be defined by 

v l (x) := v(exi, £X 2 , e(i + x 3 )). 
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By Lemma 14.31 there exists a positive constant c such that for every i, 

I |Vu ®| 2 dx < c( I \Ev l \ 2 dx + f |Vi/| 2 dxV 
Jm ' Jm Jf ' 

and changing variables, we find 

f | Vz >| 2 dx < c( f \Ev\ 2 dx + I |Vf)| 2 da:V 
Jm* ' Jm* Jf* ' 

where the constant c does not depend neither on i nor on e. Summing over i we 
find 

I | V4>| 2 dx < c( [ | Ev | 2 dx + f |Vf)| 2 dx'), 

J M e ' J M e Jf s ' 

and changing variables, using again the notation introduced in Remark 12.11 we 
conclude the proof. □ 

Lemma 4.5. There exists a constant c > 0 such that for every v £ H dd (F) we 
have 

Wv\\ L ^ F) + ||( j R £ )- 1 u|| l2(f) < - £ \\E^v\\ L 2 (F) . 

Proof. This follows immediately from a rescaled Korn’s inequality proved in [2|- 
Indeed, using the notation introduced in Remark [2.11 let v = ( R e )~ 1 v o (r E ) _1 . 
Then in [2] it has been proved that 


[ |Vf)| 2 + |D| 2 dx < f \Ev\ 2 dx. 
d Fp ^ J F e 


>F e JF e 

The result follows by simply changing variables. 

Lemma 4.6. For every v £ H dd (Q) we have 

F e (v) > c\\£H e v\\ 2 L 2 (M) . 

Proof. Since C is uniformly positive definite we have 

F e (v) > m(\\E £ v\\ 2 L 2( F ) + £ 2 ||-E e x || 2 2( M >) 

> c£ 2 (||7L‘"u|||2( F ) + ||-E‘"HI|2(m))i 

> c£ 2 \\H e v\\ 2 L2{M) 

where we first used Lemma [43] and then Lemma [4T4l 


□ 


□ 


Lemma 4.7. For every v £ H dd (Q) we have 

^M> C (||u'|| 2 ffl(M) + ||u3|| 2 2(/;ffl(ajViS)) ), 

with v 1 = (vi,V 2 )- 

Proof. From Lemma H31 Q, Poincare inequality, and small e, we have that 

JF e(v) > c([|— + Wfov'Wl^M)) — c ll^ 7 l;/ |li 2 (M) — C \W lljH 1 (M) • 
Also, from Lemma 14.61 and (foil we find that 
(24) F E (v)>c\\Vv 3 \\ 2 LHM y 
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From the following Poincare’s type inequality 


|«31 2 dx < c 


( / |V'a 3 | 


dx 


J m Jn 

Lemma 14721 and (1241) . it follows that 

T e {v) > c||u 3 ||| 2(M) , 

and hence the statement of the Lemma follows. 



□ 


Lemma 4.8. Let u e £ H dd (fl) be the solution of problem Q. Then 

sup F e (u e ) < +oo. 

£ 

Proof. By taking <p = u £ in problem o, we find 

Fe{u e ) < ||/||L2(n)l|w £ ||L2(n) 5 
while from Lemma 14.21 and Lemma 14.71 we deduce that 

C IKIIis(fi) < Fe(u e ). 

The Lemma follows from the previous two inequalities. □ 


From Lemmas rot m rra rrn and Ol one deduces (HDD , d). and (1221) : thus, 
the proof of Theorem 14.11 is achieved. 

We now prove some compactness properties using the previous a priori estimates. 

Lemma 4.9. There exist ( u , v,w) SWxVxW and z 3 £ L 2 (I\ H l (u})) with z 3 = 0 
almost everywhere in F, such that the sequence of solutions u e of problem 0 fullfills 
(up to a subsequence) 

u e a u a in H 1 (LI), 
u 3 u 3 + z 3 * n L 2 (7; JT 1 (cv)). 

Moreover, 

( E e u e ) 33 XF ( Eu) 33 xf, £{E e u e ) 33 XM 0, e(E e u e ) a3 XM ^ d a z 3 x M , 

( E e u e ) a3 XF ( Ev)a3XF , ( E e u e ) a/ 3XF (Ew) a pXF in L 2 (fl), 

in L 2 (Ll). 


Proof. From (1221) . up to subsequences, it follows that 

u% u a in H 1 ^), u% u 3 in L 2 (I; 

for some u a £ H dd (Q) and u 3 £ L 2 (I- 1 H 1 (uj)). From (0 and (1221) we deduce that 
S7'u e a —X 0 in L 2 (fi;K 2 ) thus Vu a = 0. Hence there exist £ 4Lq( 0,£) such that 


u a (x) = f, a {x 3 ) for a.e. x £ LI. 

From 4 HD it follows that u 3 —*■ u 3 in H 1 (F) and since ( Eu e )i a —X 0 in L 2 (F), we 
have that (u\,U 2 , u 3 ) £ BNdd{F), that is: £ H 3 (0, £) and there exists a function 

( 3 £ Eq(0,£) such that 


u a =£ a , = $3 -x a ff a in F. 

Set 

u 3 := f 3 - x a £' a £ H dd (Ll), z 3 := u 3 - u 3 £ L 2 (I; H 1 ^)). 

Then, z 3 = 0 in F, u := (iti,U 2 ,u 3 ) £ BNddiLl), and u 3 u 3 + z 3 in L 2 (I; H 1 ^)). 
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Since (E e u e ) 33 = d 3 u% we have that 

(■ E e U e )33XF —" 8 3 U 3 XF = 8 3 U 3 XF = {Eu) 33 XF 

in L 2 (fl). Since e(E e u E ) 33 = e(H £ u E ) 33 from (I2T1) it follows that e(E £ u £ ) 33 is 
weakly convergent in 1?(M ), but since eu\ —>• 0 in L 2 (M) we deduce that 

e(E e u e ) 3 3XM ^ 0 in L 2 (Q). 

Finally, 

e{E e u e )a3XM = -^{8aul+d3U e a )xM -^(8 a u 3 +d a z 3 +d 3 u a )xM = ( Eu)a3+-8 a Z3 

in L 2 (Q). The stated convergence result follows by recalling that u £ BN dd(f2) and 
hence that ( Eu) a3 = 0. 

For the proof of the last two convergences, see [15] and □ 


In Lemmas 3.1 and 3.2 of m it is shown that 
(25) e(H e u e ) 2 iXF Vxf in L 2 {M), 

where ft £ Hq(0,£)-, moreover, with v as in Lemma HUl we have (Ev) a3 = ^(ft'x^ + 
d a v 3 ). 

Lemma 4.10. Let ft £ Hq(I) be as in (|25j) . There exists z a £ L 2 (/;lL 1 (w)) with 
z a = 0 almost everywhere in F, such that up to a subsequence, 


£ 


- -I — dx ftx„ + z a in L 2 (/; H 1 (w)). 

JUJ £ 


Moreover, 


Proof. Let 


e(E E u E ) a/ 3XM -*■ ( Ez) a pXM in L 2 (Ll). 


Z « := T~ 


/. 


-2-dx' 

£ 


Since + z E a dx' = 0, by Poincare’s inequality we have 

J UJ 

[ \z E a \ 2 dx'<C f \Vz e a \ 2 dx' = c\ I \V'u E a \ 2 dx'<C f \£H E u e \ 2 dx', 

J UJ J UJ £ J U} j UJ 

a.e. in I = (0,£). Integrating in x 3 over I we deduce, thanks to (]22j) . that z E is 
bounded in L 2 (I; H 1 (w)) so that up to a subsequence, 

z a in L 2 (/;IL 1 (w)), 

for some z a £ L 2 (/;IL 1 (w)) such that T z a dx' = 0. From (EUl) we have that 

J UJ 

eE e u E —> 0 in L 2 (F-,M. 3x3 ) and hence, taking into account (l?5l) . we have that 

e{H E u E U^(^ ° ) inL 2 (F;R 2x2 ). 

But e(H £ u £ ) a p = -8pu E a = dpz E and hence 

d\Zi = d 2 z 2 = 0 , diz 2 = -d 2 z\ = ft in F. 

By integration we find 

Z a = ftXa + d a , in F, 
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with d a £ E 2 (I). But since j- z a dx' = 0, we deduce, by invoking ©, that d a = 0. 

J UJ 

We set 

z a '■= z a — infb 

Then z a = 0 in F and £ L 2 (/;i/ 1 (w)) . Also, 

£{E e U e )aPXM = —(dallp + d/3U e a )XM = ( Ez E ) a 0XM - 1 ‘ ( Ez)aPXM = {Ez) a pXM 


L 2 (n). 


□ 


Proof of Theorem 13.11 The convergences stated in Theorem 13.11 with strong 
convergence replaced by weak convergence are proved in Lemmas 14.91 and 14.101 At 
the end of the proof we will show that indeed are strong. To identify the limit 
problem take test function p in © of the following form 

Pa (x) = U a ( x 3 ) + £V a (x) + £Z a (x) + £ 2 U) a (x) 

¥> 3 (x) = u 3 {x) + Z 3 (x) + £V 3 (x), 


(26) 

with 


u = ( u a ,u 3 ) £ BAfdd(tt), v = {v a ,v 3 ) £ Tlddiyi) niL] d (fl) 

and 

w = (w a , 0) £ 7 ZV^iF) n H dd (n), z = (z a , z 3 ) £ H dd (n), with z t = 0 in F. 
Then 

i^E (p^atf3 — T (Ew^cxfi : 

( E e p) a3 = — d a Z 3 XM + \d 3 ZaXM + (Ev) a 3 + £{Ew) a 3 , 

Ze Z 

{E E p) 33 = (Eu ) 3 3 + ( Ez) 33 Xm + £{Ev) a3 - 

Taking p in © as in (l26|) . letting £ go to zero, and using the density of 77 dd (f2) fl 
iL dd (fl), 7 ZT> 2 (F) (~l H dd (fl) and Z fl H dd (fl) into V, W and Z respectively, we 
deduce (fl3l) . 

We now prove the strong convergences stated in the theorem. Set 

_ f (Ez'jafi ~2^aZ 3 

\daZ 3 0 


Ef ■■= 


( {Eu>)ap ( Ev) a3 


Em, := 


(Ev)a 3 (Eu) 33 
Taking p = u e in © and passing to the limit we find 

lim [ f xf + £ 2 Xm)CE e u e ■ E e u e dx = lim [ f ■ u e dx 

E_ >° Jn e_> ' 0 Jn 


f a U a + f 3 {u 3 + Z 3 )dx 


J 

Jn 

/ CEf ■ EfXF + <EE m ■ E m XM dx. 
Jn 


Since the norm A K > f n CA ■ A dx is equivalent to A || Al|| z. 2 (o) ^ it follows that 
J™ IKxf + £XM)E E u e \\ L 2(Q,) = \\xfE/ + XME m \\L 2 (n), 

which, together with the weak convergence, implies that 


(27) 


(xf + £Xm)E e u e ->• XFEf + XMEm in L 2 (tt). 
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We now prove the strong convergence of the displacement components. By Korn 
inequality we have, for e small, that \\E £ u £ \\l 2 (f) > \\E uE \\l 2 (f) > c||ti e ||#i(f) 
and hence, by means of (HZD, we deduce that u £ is a Cauchy sequence in H 1 (F). 
By means of Lemma ( 14 . 5 J) first and then Lemma ( 14 . 41 ) we deduce that eH £ u £ is 
a Cauchy sequence in L 2 (fl) and then, with the same argument used in the proof 
of Lemma 63), we find that u £ a is a Cauchy sequence in H 1 (M) and that itf is 
a Cauchy sequence in L 2 (/; 7 L 1 (w \ w)). This completes the proof since we have 
already shown that u £ is a Cauchy sequence in H 1 (F). □ 


Proof of Theorem 13.31 All the convergences stated in Theorem 13.31 may be 
obtained choosing a test function ip £ L 2 {Tl) that depends only on the X3 variable 
in the corresponding convergences stated in Theorem 13 .II To prove that U can be 
written as 

U(x3 ) = + 113 dx' + mo(x3)+ o /3 dx' + 77100(23), 

J UJ J ui\uj 

we first recall that, by definition, U(x 3) = -j- 113 + ^3 dx', where 23 is the solution 

J UJ 

of problem ( 1131 ) with u = v = w = 0 , that is 
23 £ Z, 


( 28 ) 


/ C(z) ( ( f^ 3)a/3 2 

J fi V l d « z 3 


1 d a z 3 

0 


^0,23 0 


XM dx = 

' f3(x)z 3 (x)dx, 


. V23 £ Z. 

Since problem (E 51 ) is linear we can write 23 as 23 = 2 1 + 2 00 , where 2 00 is the 

solution of problem (E 51 ) but with fs(x) replaced by fs(x) — + fs{x', X3) dx', and 

J u)\uj 

2 1 is the solution of problem (l 28 l) but with fs(x) replaced by j- fs(x',x 3) dx'. 

J uj\uj 

By linearity and uniqueness one can check that 2 1 = j- fs(x',x 3) dx'z°, where z° 

J uj\lu 

is the solution of problem (TT 4 l) . The representation of U then immediately follows. 
□ 


5 . The homogenization problem: proofs of the main results 

We start by proving some a priori estimates. 

Lemma 5.1. The sequence of solutions u £ of problem 61 satisfies the following 
a priori estimates 

( 29 ) sup||u||| L 2 (n) <+00, 

£ 

( 30 ) sup||£u £ || ff i (n) <+00. 

£ 

Moreover, setting 

J e (u e ) := f {xF e + £ 2 XM e )\Eu e \ 2 dx, 

Jn 
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we also have 

(31) sup J e (u e ) < +oo. 

£ 

Proof. By Korn inequality, there exists a constant C such that 

f I |u| 2 dydx 3 <C I I \Ev \ 2 dydx 3 £ H\ d {D x I). 
J I J D J I J D 


With v as 

v a {y,x 3 ) = eu E a {ey + ei,x 3 )\ v 3 (y, x 3 ) = u E 3 (ey + ei, x 3 ), VyeDxI, 

we obtain, after the change of variable x := (ey + ei, x 3 ) £ D\ x I := FI, that 

[ (s 2 \u %\ 2 + \ul\ 2 )dx <C I ^{e^Eu^ 2 + e 2 \{Eu E ) a3 \ 2 ) + \(Eu E ) 33 \ 2 dx 
J n JF 'e a p 

< C f \Eu e \ 2 dx. 

■>n 

Taking the sum over i in the previous inequality, we get the estimate 


(32) 


[ \u 3 \ 2 dx < C f \Eu e \ 2 dx. 
JF r J F P 


We now apply the following Poincare’s type inequality 

|V'w| 2 <iE+ f \v\ 2 dx), 


Iwl 2 dx < 


IM 


with fl := 7 x J, F := D x 7, M := (Y\D) x I and v(y, x 3 ) = u§(ey + ei, x 3 ) so 
that using the same change of variable as above and the classical Korn inequality 
||£M e ||#i(f 2 ) < C\\eEu E \\L 2 (Q), we infer, after taking the sum over i, that 


(33) 


f \u 3 \ 2 dx < C f \eEu e \ 2 + j |m || 2 dx. 
Jm, Jq Jf, 


Since by definition J e {u e ) = / (\F e F e 2 XM e )\Eu E \ 2 dx, we deduce from (l32l) and 

Jn 

(l33l) that 


(34) 


[ |m || 2 dx < CJ £ (u e ). 
hi 


We now take ip = u e in (1121) . Applying the Cauchy-Schwarz inequality in the right 
hand side and using the strong ellipticity of C together with (|34D and the Korn 
inequality ||£w e || ff i ( n) < C||£.Eu E || L 2 ( n ), we get 

(35) _ 

Je{u E ) < C II eu s a || L 2 (n) Tc || u% ||L2 (fi) < c II eu e a || H i(n) +c || u% ||z, 2 (f2) < c^/j e {u e ). 

From (l35l) . we first deduce (1311) and then, by means of (l35l) . we deduce (l29l) and 

□ 
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Lemma 5.2. There exist (ui,u) l ,'d) £ U h , v 3 £ V h , and w £ W h such that the 
sequence of solutions u e of problem da fullfills (up to a subsequence) 

(36) eu e a ->■ u a in 

(37) eV'< — Vu a + Vy a , 

(38) u% u 3 , 

(39) eV'u§ V' y u 3 . 

Moreover, setting v a := y^d, we have 

(■ E«w) a p [Ey*v) a 3 

( E «* v) a 3 (£11)33 


(£u e )xf s 

e{Eu e )xM e 


(Eu + Eyu 1 )^ {Ey*u) a 3 
(£ yx w)a3 0 

where the operator E yx has been defined in (ED- 


Xd{v), 
Xy\d{v), 


Proof. By ([29]) . (f3fl)l . and classical compactness results related to two-scale con¬ 
vergence, recalled in Section [2j there exist u a £ H dd (Q), u\,u 3 £ L 2 (fi; H#(Y)) 
such that, up to subsequences, the convergences stated in (1361) (l39l) hold. From 
m it follows that 

(40) sup \\xF e Eu e \\ L 2 (n) < +00 sup \\£Eu e \\ L 2 {n) < +00. 

£ £ 

Thus, up to subsequences, we have that 

XF e £(Eu e ) a/ 3 0, 

while taking into account (1371) it also follows that 

XF e e(Eu e ) a p ((Eu ) at 3 + (£ y u 1 ) a , s )xD(?/)- 

Hence 

(Eu) a p + ( E y u 1 ) a p = 0 for a.e. (x,y) £ x D. 

Since u does not depend on y , from the above equation we deduce that 
ul{x,y) = -(Eu{x)) a0 yp+d{x)y^ £ L 2 (fi; H^(D)) 

for some function i9 £ L 2 (fl). We remark that for every x the sum r d(x)y(( denotes 
an infinitesimal rigid rotation. Since 

x a . d 2 ux - diu 2 R 

\Eu) a pyp = dpu a yp H--- y a , 

we can write 

ui(x,y) = -dpu a (x)yp + d(x)y(( for a.e. (x,y) £ ft x D, 

where we set 

„ a &2U1 — <9i1t2 

- - - . 

By considering the component ( Eu e ) 33 , and (l40j) we get that 
d 3 u%XF e K 33 (x,y)xD(y), 

for some function K 33 £ L 2 (fi x D ). For <j> £ 23(H) and if £ C#(Y), we have 

(41) f d 3 u e 3 XF e H x )^(—) dx = - f u 3 XF e d 3 <f>(x)i/>( —) dx 

Jn £ Jn £ 
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and by passing to the limit we find 

/ / K 33 {x,y)cp{x)ip{y) dxdy =- / u 3 (x,y)d 3 <p(x)ip(y) dxdy. 

J Cl J D J Cl J D 


Therefore, u 3 £ L 2 (u> x D\H l (I)) and d 3 u 3 = K 33 a.e. in fl x D. In addition, we 
have u 3 £ L 2 (u> x D\Hq(I)). Indeed, identity (14T1) holds also for test functions 
cp £ C°°(fl) and ip £ C#(Y), since £ H dd (fl), and hence we deduce that 


/ / d 3 u 3 XD(y)(p(x)ip(y) dxdy 


U3{x,y)XD{y)d 3 (j>(x)ip(y) dxdy. 


JQJY Jujy 

Integrating by parts the integral on the right, we conclude that 

(42) u 3 (x',0,y) = u 3 (x',i,y) = 0. 

From (HD1) we deduce that 


2e(Eu E ) a3 XF e 0, 
while, by means of (1551) and (1351) we conclude that 

2 e{Eu e ) a3 XF e = (ed a u e 3 + ed 3 u e a )xF e (d v a u 3 + d 3 u a )XD{y)- 

In the above limit we used the fact that ed 3 u E a —“■ d 3 u a since there is no fast 
variable in the X3-direction. Thus 


(d> 3 + d 3 u a )xD(y) =0, 

from which we easily deduce that 


u 3 (x, y) = -y a d 3 u a (x) + u 3 (x), for a.e. (x, y) £ II x D, 

for some u 3 £ L 2 (uj ; Hq(I)). Note that this equality implies that u a also belongs to 
the space L 2 (w; Hq (/)). 

We now prove the convergent of the rescaled strains. Thanks to (l3lT) . there exist 
symmetric matrix fields K,G £ L 2 (fI;R 3x3 ) such that, up to subsequences, 


(43) Eu e XF e —^ K X D{y), 

(44) eEu e XM € Gxy\d{v )■ 

We now characterize the components of K and Z. That ^33 = d 3 u 3 a.e. in 
x D has been already proven above. The proof of the characterization of the 
two-scale limits K a 3 = ^(y^;d 3 d + d y v 3 ) and K a p = (E y w) af 3 for some v 3 £ V h , 
1 ? £ L 2 (uj] and w £ W h is very similar to the one given in [5j, Propositions 

4.1 and 4.2. For the sake of brevity, we will refrain in reproducing it here. 

We now prove that d 3 d = d 3 d in the sense of distributions, where d is the 
function entering in the definition of zt 3 in id x D. 

For every tp £ V(VL x Y) we have that 


r „/ /* / / 

/ e(d 1 u E 2 -d 2 u\)d 3 y(x, — )dx = 2e / (Eu E ) 23 di(ip(x,—))-{Eu £ )i 3 d 2 ((p(x,—))dx. 
Jn £ Jn £ £ 


Taking p(x,y) = ip(x)<p(y) with ip £ V(p) and (p £ V(D) and evaluating the 
derivatives on the right-hand side of the identity above, by means of (l37l) and (l43l) . 
we may pass to the limit to find 


IflxD 


(diu 2 + d\u\ — d 2 u\ — d 2 u\)(pd 3 ip dxdy = 2 / (K 23 dfcp — K\ 3 d 2 cp)ip dxdy. 


IQxD 
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Taking into account the structure of u} a and of K a 3 , the identity above reduces to 


/ 2 d 4 > dxdy = 2 / — 831 9 <p ip dxdy, 

/QxD JftxD 


that is: 83d = 831) in the sense of distributions. Hence 1 ? belongs to L 2 (w; H 1 ^)) 
and not just to L 2 (fi). This implies that £ L 2 (fi; H#(Y)) fl L 2 (fi; H^D)) H 
L 2 (u} x D\H l (I)). We can assume to be equal to 0 at X 3 = 0, since, as it 
follows from its definition (EH. ir a is defined up to a function of X3. To remove 
this ambiguity we require u * . This forces d = 0 at X 3 = 0, and as a consequence 
we deduce that d = d € L 2 (w; Hq (/)). This, in turn implies that also u l a £ 
L 2 {u x D] Hq(I)). 

Thus, to sum up, we have 


(45) 


K = 


( E«w) a p 
l(y^ + dy a V3) 


l(y^ + dV a V3) 
831X3 


We now analyze (l44l) . The component G 33 is equal to zero since s{Eu e ) 33 XM e — 
—*• 0 by (|38l) . For the component G a 3 arising in ([44]) . we study the limit of each 
term in the sum 

2 e{Eu e ) a 3 \M e = (ed a u e 3 + ed 3 u e a )xM e - 
For the second term we use (l3f>l) (which implies that eu e a —^ u a ) to get 


(46) 


f x f f x f 

/ Ed 3 U e a <t>{x)il)(—)xM e dx = - £u e a d 3 (t> ip(—) xm s dx 

Jn £ Jn e 

/ / u a d 3 (j)(x) ip(y)xY\D dxdy 
Jn Jy 


—> - 


/ / 83 U a (j>(x) ijx(y)xY\D dxdy , 


Jn jy 

for all (f> £ £>(f2) and ip £ C#(Y). By virtue of (13^1) . the two-scale limit of 580 . 1 X 3 is 
equal to 8 ^ 1 x 3 , thus 


2e(Eu e )a3XM e —^ {da u 3 + &3U a )XY\D = G a 3XY\D- 
To identify the components G a p we use Eil) to get 

2e{Eu e ) a pXM e = (edaUp + edpu%)xM e 

(d a up + 8 v a up + 8pu a + d^ul)xY\D 
= 2 (Eu + E v u 1 ) a pXY\D = 2G a pXY\D- 


Hence 



and the Lemma is proven. 


□ 


Proof of Theorem 13.51 The convergences stated in the Theorem, with the strong 
two-scale convergences replaced by the weak two scale convergence, have been al¬ 
ready proven in Lemma 15.21 Set 

cr e := C Eu e , t s := eCEu E 
then, by ED and (l3ll) we deduce that 

&Xd, t s xm s txy\d, 


a Xf, 
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for some <Jij,Tij £ L 2 (fl x Y). For later use, we now prove that 


(48) 


ID 


v a p(E v ip) a p + a a 3 dli/j 3 dy = 0 


almost everywhere in fl for every ip £ H^{Y). Indeed, take p(x) = cp(x)ip(x'/e), 
with (p £ V(fl) and ip £ H#(Y). in (fl2l) to find 

f i f x f x f 

/ (& E XF e + £ T e XM e )-(tpE(p+(p-E y ip)dx = / ef a (x,—)-ip a (x)+f 3 (x,—)-ip 3 (x)dx. 

Jn £ Jn £ £ 

Multiplying the above identity by e and by passing to the limit we deduce 


[ [ <r(x, V) ■ {4>(x)E v ip(y))xD dydx = 0, 
Jn Jy 


which implies C51) . In particular, C51) implies that 


(49) 


/ cr a p dy= a a py 7 dy = 0 
J D J D 


as it can be checked by taking, in (H51) . ip a £ H\(Y) and equal to ip a = a < J l y 1 + 
bjgV-yds on D, where a“, b‘J s £ {0,1}. Also, from (l48l) we deduce 


(50) 


cr a3 dy= 0-13J/1 dy= a 23 y 2 dy= cr 13 y 2 + a 23 y x dy = 0. 


id 


id 


ID 


ID 


We are now in a position to exhibit the appropriate test function to pass to the 
limit in (1121) . We set 


(51) 


<Pa(x) = \u a {x) + U^X, £-) + EW a (x, |-), 
V>3(x) = u 3 (x, ^r) + £V 3 (x, ^), 


where u a £ V(fl), u^,u 3 £ V(fi; H^(Y)), and w a ,v 3 £ D(Q) x T>(D). Moreover, 
we assume that 


(52) ui{x,y) = v} a {x) - y^x) - yf}dpu a {x), 

(53) u 3 {x,y) = u 3 (x) - ypd 3 up(x), 

in SlxD, where u^,u 3 , d £ 2?(f2). With this choice of test functions we find 

£(Eip) a p = (Eu + E v v}) a p + £(Eu 1 + E v w) a p + £ 2 (E v w) a p, 

2 £{Eip) a3 = (d 3 u a +d v a u 3 ) + £(8^ +d a u 3 +d v a v 3 ) + E 2 (d 3 w a +d a v 3 ), 
£(Elfi) 33 = £d 3 u 3 + £ 2 8 3 v 3 , 

and taking into account (1521) and (l53l) we deduce that 

( Eip) a pXF e = [(Eu 1 + E v w) a p + £(E y w) a p]xF e , 

^(E ip) a 3 XF e = \(8 3 u 1 a + d a u 3 +dlv 3 )+£(d 3 w <x +d a v 3 )\xF e , 
(Eip) 33 XF e = [d 3 u 3 +£d 3 v 3 \xF e - 

Problem (fl2l) rewrites as 

r r x f x f 

/ ( a e • E(pxF e + r e ■ £EipxM e dx= £f a (x, — )■ <p a (x) + f 3 (x, —) • tp 3 (x) dx, 
Jn Jn £ £ 
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whose limit is 

<7 • 

(54) 


We claim that 




(55) 

and that 

(56) 


/ / (JapiEu 1 )^ dydx = 0 
Jn Jd 


v a 3 (d 3 u], + d a u 3 + d y v 3 ) dydx = / <r a3 (y*d 3 d + d y v 3 ) dydx. 

in Jd JnJ d 

Indeed, by means of (l52l) we compute 
{Eu l ) a p{x,y) = {Eu 1 ) a p{x) - ^(y%dp&(x) + yfd a d(x)) - y^(Eu) a p(.x) 


which, together with C^l) . implies (l55l) . To prove (15^1) . note that, 

Ua{x,y) + d a u 3 (x,y) + d y v 3 (x,y) = 

[ya d 3^( x ) + d v a v 3 {x,y )] + [9 3 u^(x) + d a u 3 {x)\ - ypd 3 (Eu) a p(x), 
where we used (l52l) and (l53l) . Taking into account (l50l) we deduce that 


/ / (r a 3 [d 3 U 1 a (x,y) + d a u 3 (x,y)+d y v 3 (x,y)\dydx 

Jn Jd 

o’a 3 [yud 3 'd{x) + d y v 3 (x, y)} - a a3 ypd 3 (Eu) a p(x) dydx 
<Xa3[ya d 3 , d{x) + d*v 3 (x, y)} - ai 3 yid 3 (Eu)n(x) 


n Jd 


n Jd 


-~cr 23 y 2 d 3 (Eu) 22 (x) - (a 13 y 2 + V23yi)d 3 {Eu) 12 (x) dydx 


and using again (15U1) we obtain (15^1) . Thanks to (l55l) and (15^1) . (1571) reduces to 



which holds for every u a £ X>(S1), u^,u 3 £ P(f l',H#(Y)), w a ,v 3 £ £>(f 1) x V(D) : 
and d £ 2?(fi) that satisfy (152)) and (l53l) . By density (1571) holds also for (ui,u]y) £ 
U h , (—y 2 d,yid,v 3 ) £ V h , and w £ W h . 

To conclude the proof it suffices to show that 


a\D = C 


{E y w) a p (E yx v) a3 
(E»*v) a 3 (Eu) 33 


XD 
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TXY\D = C 


Xy\d- 


and 

(Eu + Eyu 1 )^ (Ey*u) a3 

{Ey*u ) a3 o 

But, these identities follow essentially from the definition of two-scale convergence 
and Lemma P For instance, with ip £ L 2 (Q ; C#(Y)) we have 


f x f f x f x f 

/ Pij{x)XF e v{x,—)dx= C ij ki(x,—)Eu% l (x)xFM x >—)dx 

Jn e Jn £ £ 


in £ Jn 

and by passing to the limit we find 


/ / <?(x,y)XDp(x,y)dydx= / C ijk i(x, y)G(x, y)xDp(x 1 y) dydx 
JnJY JnJY 

r 2 / 


where we used the fact that C iju<p £ This identity implies a\D = 

CKxd- 

The strong two-scale convergence of the rescaled strains can be proven as it was 
done in the rod problem. □ 

Proof of Theorem 13.61 With 

Za : = wi + ypd 3 u a - fiy*, z 3 := u 3 - (f u 3 dy - y a d 3 u a ) 

J D 

in x Y, it follows that the bilinear form of problem (fI71) defined on Q x (Y \D) 
rewrites as 


f f f (Eu + EVu 1 )^ (Ey*u ) a3 \ 

Jn Jy\d V (E yx u ) a 3 0 J 


(Eu + E v u l ) a p (Ey x u) a 3 

(Ey x u ) a3 o 


dydx 


f f ( (Ey z ) a p \d y a z 3 \ ( (Eyz) a0 \d v a z 3 \ 

= LLo C \ o j •(. % ) d 

By means of the definitions of aj, and u 3 on x D we infer that Zi = 0 in x D. 
With this new notation we have that 

U(x) := 4- u 3 dy = j- z 3 + (1 u 3 dy - y a d 3 u a ) dy = f u 3 dy + i z 3 dy. 

Jy Jy Jd Jd Jy 

The proof of theorem is now very similar to the proof of the corresponding Theorem 
13.31 in the rod problem. Since it can be achieved by following exactly the same 
steps, with the role of the variable x' played by the variable y and the spaces 2, Z$ 
replaced by Z h , Zq, we refrain from writing all the details of the proof. □ 
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